The simulations can be reproduced through the master script ``run_sims.R``, which sources ``estimator_functions.R`` and ``simulation_study_functions.R``. 

The figures and tables in the text can be reproduced as follows:

	- Figure 1 (Critical J): ``critical_J.R``
	- Figure 2 (Simulation RMSE): ``sim_text_figures.R``
	- Figure 3 (Simulation coverage): ``sim_text_figures.R``
	- Table 3 (HSIS site-ATE estimates): ``analyze_pseudo_data.R``
	- Figure 4 (HSIS self-efficiency): ``analyze_pseudo_data.R``
	- Table 4 (HSIS diagnostic statistics): ``analyze_pseudo_data.R``
	- Figure 5 (HSIS RMSE): ``hsis_rmse_figure.R``
	- Figure 6 (HSIS coverage): ``hsis_coverage_sim.R``

To reproduce the figures in the Supplement:

	- Section E (Simulation distribution illustrations)
    	- E.2.1, E.3.1, E.4.1: ``dist_illustrations.R``
	- Section F.1 (Asymptotic approximations)
		- F.1.1, F.1.2, F.1.3 (RMSE): ``approx_sup_figures.R``
    	- F.1.4 (Variance): ``approx_sup_figures.R``
    	- F.1.5 (Bias): ``approx_sup_figures.R``
	- Section F.2 (Accuracy of approximations)
    	- F.2.1, F.2.2, F.2.3: ``sim_sup_figures.R``
	- Section F.3 (Simulation RMSE)
    	- F.3.1, F.3.2 (hetero/nonlinear): ``sim_sup_figures.R``
    	- F.3.3, F.3.4, F.3.5 (normal/heavy): ``sim_sup_figures.R``
	- Section F.4 (Simulation coverage)
		- F.4.1, F.4.2, F.4.3: ``sim_sup_figures.R``


Running in serial on a standard desktop, the full simulations of the paper took about 3 days to complete. Multicore or cluster parallelization would be trivial. ``run_sims.R`` includes lines of code to run a subset of the simulation that takes roughly 2 hours to complete, and clear instructions on how to run the rest of the simulations (essentially just by changing a function argument). The HSIS-based simulation, in ``hsis_coverage_sim.R``, to examine coverage took roughly 50 minutes to complete. Other scripts run within minutes.